Asymmetric voltage noise in superconducting tunnel junctions with electromagnetic environment 
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We investigate theoretically the V-I characteristics and voltage noise of superconducting tunnel junctions 
with small critical current via the matrix-continued-fraction method. Special attention is paid to the large 
hysteresis in the V-I characteristics and to the voltage-noise anomaly observed in preliminary experiments. The 
current dependence of the voltage noise shows a strong asymmetry between the forward and backward current 
ramping and a discontinuous change of the noise close to its measured maximum occurring at the switching 
current. We show that both the large hysteresis and the voltage noise anomaly in this current-biased setup 
are a consequence of the influence of junction's electromagnetic environment. Typically, the voltage-noise 
dependence on the junction current parametrized by the external drive contains a loop which is, however, not 
observed experimentally because of the implementation of the junction-current ramp. Skipping a part of the 
loop is responsible for the observed hysteresis as well as the noise anomaly. 

PACS numbers: 72.70.+m; 74.50.+r; 85.25.cp 
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I. INTRODUCTION 

For small-area Josephson junctions, where the critical cur- 
rent is in the nanoampere range and the capacitance is in 
femtofarads, the dynamics of the Josephson phase is strongly 
influenced by the electromagnetic environment |Q]] . This can 
have profound effects on the measured V-I characteristics and 
voltage noise. For example, depending on the realization of 
the environment the same junction can be in the underdamped 
regime, with a typical hysteretic behavior of the V-I curves 
between the forward and backward source-current ramping, 
but also in the nonhysteretic overdamped regime. 

A recent preliminary experimental study by Pertti Hako- 
nen's group in Helsinki |2j[] focusing on the voltage noise 
of small-area tunnel junctions revealed a strong forward 
vs. backward current-ramping mean-voltage and voltage- 
noise asymmetry and hysteresis (see Fig. [TJ. Furthermore, a 
discontinuous change of the voltage noise coinciding with the 
switching current was observed. The purpose of this paper 
is to show that such surprising effects can be understood as a 
consequence of junction's electromagnetic environment. 

The organization of our work is as follows. First, we in- 
troduce in Sec. [TT] a simple model circuit close to those used 
in experiments for obtaining the critical current of real junc- 
tions and shortly discuss the matrix continued fraction method 
used to analyze properties of the model circuit. In Sec. fXTTJ 
we present the basic ideas of our explanation. Finally, we dis- 
cuss our numerical results and show that they are in qualitative 
agreement with the experiment in Sec. [TV] 



II. CIRCUIT AND METHOD 

As it is known the electromagnetic environment can be en- 
gineered to place any Josephson junction in the overdamped 
regime [3]. We investigate the influence of such environments 
on the junction properties in more detail. The considered 
circuit is shown in Fig. [2j where a Josephson junction (red 
square) with critical current I c , intrinsic capacitance Cj, and 



.x10 




FIG. 1 : Experimental dependence of the mean junction- voltage on 
the bias current (top panel) with the corresponding voltage-noise (in 
arbitrary units) dependence (bottom panel). Black line represents 
forward current ramping while the red corresponds to the backward 
direction. Unpublished data courtesy of P. Hakonen Q]. 



resistance Rj is biased by a circuit with a capacitor C, resistor 
R, and an ideal voltage source V s . This circuit is equivalent to 
a current source I s = V s /R in parallel with R, C, and the junc- 
tion. Identical circuit was used to show that the commonly 
observed reduction of the maximum supercurrent in an ultra- 
small junction is not an intrinsic junction property, but is due 
to its electromagnetic environment [3[]. Its advantage lies in 
the fact that, in contrast to the usual current-bias method (ap- 
plied directly to the junction) giving access only to the positive 
differential resistance part of the V-I characteristic, all points 
on the V-I characteristic are in principle achievable. That's 
because both the current / through the junction (measured by 
ideal amperemeter /) and the junction voltage V (measured by 
ideal voltmeter V) are average quantities adjusted to the global 
drive. 
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FIG. 2: Idealized circuit for measurement of the supercurrent of 
Josephson junctions. The red square contains a junction modeled 
by the RCSJ model. 



The dynamics of this circuit following from Kirchhoff's 
laws and Josephson's relations is described by the Langevin 
equations 
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where two current-noise sources £,r and £,r, associated with 
(mutually uncorrected) thermal fluctuations in resistors R and 
Rj are assumed to be simple Gaussian white noises satisfying 
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It is useful to introduce the equivalent parallel resistance 
Rp = RRj/ (R + Rj) and the equivalent parallel capacitance 
C p = C + Cj. Then the quality factor of the circuit reads 



Q = (OpRpCp, 



(3) 



where co p = y/I c /q>oC p is the plasma frequency with (po = 
h/2e the reduced flux quantum. For the purpose of theoreti- 
cal analysis it is useful to rewrite Eq. ([TJ into the dimension- 
less form by introducing the following dimensionless quan- 
tities |01 - junction voltage v = YW' J unct i° n current ; = j-, 
time T = &pt, temperature ® = (with £g the Boltzmann 
constant), driving force i s = jf-, and, finally, the composite 
Gaussian white noise £ with the correlation function 

<C(xi);(x 2 )) = ^<(^i)-^(*i))(^(fe)-^(/a))> 

= 2708(1! -x 2 ), (4) 

where y = 1 /Q is the the damping coefficient. 

Using these definitions Eq. (Q]i can be reformulated in the 
dimensionless form (v(t) = 3(p/dx) 



dv(x) 
3x 



= i s -yv(x) - sincp(T) + £(x). 



(5) 



This equation is formally identical to the RCSJ model fH, 
which is commonly used for description of Josephson junc- 
tions without the environment and, thus, the methods devel- 
oped for the RCSJ model can be applied. We have used the 
matrix continued fraction method Jg, Sec. 11.5] for the solu- 
tion of the Fokker-Planck equation associated to the Langevin 
Eq. @ for the probability distribution function W(<p, v,x) 
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The probability distribution function is obtained numerically 
by expanding the v-part of W((p, v,x) into quantum oscillator 
basis functions, thus obtaining a tridiagonal coupled system 
of differential equation, and the 27t-periodic (p-part into the 
Fourier series Ja, Sec. 11.5]. We restrict our analysis here to 
the stationary case where the average value of junction voltage 
v and junction current i can be computed as 
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(v) = J d<p J dvvW(<p,v,°o), 
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(7) 
(8) 



where the environmental parameter p =R/Rj, and the voltage 
autocorrelation function reads Sec. 7.2] 



In 



(v(x)v(O)) = j dcp J dvve" |T l Lj7i W((p,v,°o) 



(9) 



Finally, the steady state voltage noise is obtained by time in- 
tegral of the (connected) voltage autocorrelation function 



S = 



J dx((v(x)v(0))-(v(x))(v(0))). 



(10) 



III. BASIC MECHANISM 

In Fig. |3la) we plot a typical curve of the mean junction 
current (i) dependence on the mean junction voltage (v) for 
temperature = 0.04, quality factor Q = 1, and different en- 
vironment parameters p. Here, limit p — > °° corresponds to the 
case without the environment. For that case the mean current 
(i) is a monotonic function of the voltage (v). This changes 
with the introduction of the environment when a local maxi- 
mum and local minimum (decreasing with decreasing p) of (z) 
develop. Fig. |3b) shows the voltage dependence of the volt- 
age noise — it turns out that for a given Q this quantity does 
not depend on the environment details and, thus, all curves for 
various p coincide. 

The non-monotonic dependence of current on the mean 
voltage has serious consequences for current-biased junctions 
as illustrated in Fig. St a). The mean voltage is not a function 
of the junction current since in the range between the current 
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FIG. 3: Dependence of the junction current (;') (a) and the voltage 
noise 5 (b) (all curves overlap) on the junction voltage (v) for = 
0.04, Q = 1 and different p's. (c) Corresponding dependence of the 
voltage noise S on the junction current (i) contains for finite p a loop. 
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FIG. 4: (a) Sketch of a typical current (i) vs. voltage (v) character- 
istic for the RCSJ model of the circuit in Fig. [2] with shown hystere- 
sis between forward and backward current ramping directions, (b) 
Corresponding picture of a typical voltage noise dependence on the 
average junction current (i) with the hysteresis and discontinuities. 



minimum and maximum 3 different voltage values can be at- 
tributed to the same current (green line in Fig.[4|a)). In prac- 
tice, when the junction current is ramped up, it will follow the 
stable branch until it hits the point A when a step back would 
be required to keep on the average characteristic (grey curve). 
Instead, the mean voltage will switch to its other stable value 
at point B (red dots). For the reverse ramp, the analogous sit- 
uation happens between points C and D (black dotted curve). 
Depending on the distance between the current extrema the 
hysteresis can be large. 

The fact that (v) is not a function of (/) has also a profound 
effect on the voltage noise dependence on (i). Typically, this 
dependence (understood as a parametric curve of both quan- 
tities parametrized by the mean voltage) contains a loop as 
shown in Fig. [3c) and further illustrated in detail in Fig. |4jb) 
by the solid grey line. There the points A, B, C, D corre- 
spond to the points in Fig. SJa). For the same reasons caus- 
ing the hysteresis in v-i characteristics one will not observe 



the whole loop in the noise either under the junction-current 
ramp conditions. The measured noise will leave the loop in 
the moment when a step back in the current bias (negative 
differential resistance) is needed to remain on the loop. This 
coincides with the points A and C, respectively, i.e. with the 
switching currents for opposite ramp directions. Therefore, 
we suggest this mechanism as responsible for the discontin- 
uous changes of the noise observed in the experiment close 
to the switching currents (see bottom panel of Fig. [TJ. Fur- 
thermore, the noise loop is not symmetric and, therefore, the 
height and the position of the peaks depends on the direction 
of the measurement as illustrated again by the black and red 
dotted lines (same color-code as in Fig. [TJ. As we will show 
later, for some parameters the lower peak becomes negligible 
and in this case we get qualitatively the same noise behavior 
as in the experiment. To sum up, we have strong evidence that 
both the large hysteresis and the asymmetric voltage noise are 
a consequence of the influence of electromagnetic environ- 
ment on the current-biased junction. 



IV. RESULTS AND DISCUSSION 

In this section we show that also other properties of the V-I 
curves and noise are qualitatively in agreement with the exper- 
iment. The junction used in the experiment has Rj = l.lkQ., 
I c = 19.8nA, and Cj = 5fF and the environment was modeled 
by R = O.lkD. and C = 0. 1 pF in all presented examples. 

In Fig. [5] we plot the V-I characteristics for different tem- 
peratures and the current dependence of the voltage noise for 
the same parameters is plotted in Fig. [6] The voltage noise 
is expressed in units of the thermal Johnson-Nyquist noise 
St = IksTRp of the equivalent parallel resistor R p . In all 
figures the solid grey line represents the full static solution 
of the circuit Fig. [2] and dotted red and black lines show the 
supposed measured curves when the junction current ramps 
are considered. From the V-I characteristics one can see 
that the measured hysteresis can by rather large and increases 
with the decreasing temperature. For the chosen parameters 
the upper drop in the mean junction voltage is approximately 
O.lmV '. It is of the same order as in the experiment (Fig. [TJ). 
There is one difference, however. While the asymptotic V- 
I lines are almost horizontal in the experiment (Fig. [JJ they 
follow the V = Rjl line in our theory. We believe that this 
discrepancy would be fixed by using a nonlinear model for 
the junction. The nonlinear model gives a better descrip- 
tion of real junctions for the voltages close to the gap voltage 
Vg 10 Sec. 2.3.2]. A simple piecewise-linear approximation 
(Rj(V) = R (a) at |V| < V G , Rj(V) = % at \V\ > V G with 
Rfa) ^ R(b)) should be enough to show that these measured 
horizontal top and bottom V-I lines have their origin in the 
almost infinite slope change of / close to the V = Vg known 
from typical V-I characteristics of a Josephson junction [7, 
Sec. 2.3.2] (note that a vertical line in the I-V curve is hori- 
zontal in the V-I characteristic). 

Our calculated voltage noise properties are also in good 
agreement with the experiment 0] (experimental data are un- 
published yet). With decreasing temperature the location of 
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FIG. 5: Junction voltage V vs. junction current / for varying temper- 
atures (for other parameter values see the main text). 



demonstrated in Fig.|7]we are not dealing exclusively with the 
thermal noise associated with the combination of the junction 
differential resistance Ro = dV jdl and resistance R. Fig. [7] 
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FIG. 6: Analogue of Fig. [5] for the voltage noise S normalized by 
S T =2k B TR p . 



the higher asymmetric noise peak shifts farther away from the 
zero current and the peak becomes steeper. The lower peak is 
very small even for low temperatures. The very same trends 
are observed also in the experiment 0]. A closer look at the 
experimental data in Fig. Q] reveals also a small drop in the 
noise close to (/) = (marked with the green circle). The re- 
alistic environment present in the experiment was apparently 
more complex than the one assumed here and, thus, quantita- 
tive differences are expected, yet our qualitative explanation 
should be valid also for more realistic environments. 

An important question is the origin of the voltage noise. As 



FIG. 7: Voltage noise 5 normalized by the equivalent noise So = 
2ksTRo of junction's differential resistance vs. the mean junction 
current / for wide range of temperatures (for other parameter values 
see the main text). Normalized-noise unitary value at zero junction 
current reflects the fluctuation-dissipation theorem. 



shows the voltage noise S divided by the So = ^ksTRop, 
where R Dp = RR D /(R + R D ). For (/) = the ratio S/S D ap- 
proaches one (as expected from dissipation-fluctuation theo- 
rem) but it exceeds one in the region where the hysteresis of 
V-I curves is observed. The origin of this extra noise is in 
the double-valued nature of the voltage (dichotomous switch- 
ing) in this regime as can be deduced from the previous RCSJ 
model studies 18[] as well as from the fact that the maximum 
of the noise is decreasing with increasing temperature (typical 
for dichotomous noise). 

In conclusion, there are several possible methods for ver- 
ification of our hypothesis. Depending on the details of the 
measurement the true nature of the V-I characteristic could 
be revealed by using the voltage ramp instead of the junction- 
current ramp. The other possibility is to use a larger number 
of measurements and slower current ramping. 
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